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1 Introduction 

The main scope of this short contribution is to call the attention on the impact 
of Negative Binomial (NB) regularity in our field as documented by the collec- 
tion of Proceedings of International Symposia on Multiparticle Dynamics since 
1972. It is also quite instructive to follow the evolution of the interpretation of 
NB regularity from its first appearance in cosmic ray physics in the sixties and 
in the accelerators region in the seventies, up to its discovery in the eighties 
in limited regions of phase space in the CERN pp energy domain, and then 
in the nineties in e + e~ annihilation at LEP in samples of events of different 
topologies up to the very recent suggestion presented at this Symposium to 
associate NB Multiplicity Distribution (MD) to quark jets of different flavour. 
This perspective is motivated by two facts. First, NB regularity in e + e~ an- 
nihilation seems to be better satisfied the more fundamental and elementary 
is the chosen level of investigation of multiparticle final states: NB regular- 
ity is less approximate for the MD of the most primitive object occurring in 
Multiparticle Dynamics, i.e., the single jet of fixed flavour, and becomes more 
and more approximate in going from single jet to sample of two jet events and 
finally to the full sample of events. Second, it should be noticed that observed 
approximations of the regularity at a certain level of investigation can be in 
part removed by properly weighting NB MD's occurring with an higher de- 
gree of accuracy in the elementary substructures present at a more primitive 
level. It is clear that if these facts would be confirmed also in other classes 
of collisions it would be quite natural to consider the single jet structure as 
an effective building block in multiparticle production as it is now in e + e~ 
annihilation, and the study of more complex final states structures as the re- 
sult of compositions of single jets whose evolution is controlled by NB MD's 
and related clan structure properties. These remarks give a deep and striking 
physical meaning to the 1979 discovery that QCD jets can be described in 
LLA in terms of Markov branching processes of NB type dominated by gluon 



self interaction, the non-linear QCD vertex, and to the subsequent observa- 
tion that clans at parton level are bremsstrahlung gluon jets. Accordingly, the 
statement that QCD parton shower formation is the main dynamical mecha- 
nism responsible for multiparticle production turns out to be a consequence 
of the intimate fractal structure of jets proposed in our field - once again - 
since 1979. Along this line of thought an integrated description of multiparti- 
cle correlations and MD's is possible. In addition generalized hadron parton 
duality as formulated in the NB framework could be indeed a highly simplified 
approximated hadronization prescription and opens a new interesting research 
line on hadronization models. 

2 Facts 

2.1 Prehistory (1966) 

NB regularity is discovered in cosmic ray physicfl. NB MD is called Polya- 
Eggenberger distribution by the name of the scientists who applied it to the 
diffusion of contagious diseases in 1926. Fluctuations of the pionization com- 
ponent in hadron showers originated by a primary hadron at different primary 
energies E are well described by NB MD's. 

Violations of multiperipheral model predictions are clearly visible starting 
from Eq ~ 30 GeV. Interestingly the Authors suggest the analytical structure 
of standard NB parameters in terms of the primary energy Eq, i.e., the average 
charged particle multiplicity n and the parameter k, which is linked to the 
dispersion D by D 2 — n + n 2 /k, turn out to be respectively 

n(E ) ~ 1.8(£ ) 1/4 (1) 

and 

k{E )~ 0.4(1 „ e ~i-8 io- 4 £o) (2) 

It follows that at E = 30 GeV, n ~ 4 and k ~ 465, whereas at E = 10 4 GeV 
n ~ 18 and k ~ 3, i.e., k is a decreasing function of the primary energy as n 
increases. A large part of the story of NB regularity in full phase space is in a 
certain sense already contained in the just mentioned results. Very large k val- 
ues suggest that the MD can be safely approximated by multiperipheral model 
predictions (a Poisson MD) and the onset of correlations among produced par- 
ticles is controlled by the increasing of 1/k at very large primary energies. See 
also the proceedings of the Wuhan 1991 symposium for an extensive discussion 
on the subject. 



2.2 NB regularity in the accelerator region (Zakopane 1972, Pavia 1973) 

NB (Polya-Eggenberger) MD is suggested for describing n-particle MD in full 
phase space in a simple model of multiparticle production which was motivated 
by the attempt to describe observed deviations from multiperipheral model pre- 
dictions in the accelerator region! NB MD is discussed at the Zakopane Sym- 
posium in June 1972. Among other interesting features it is pointed out that 
the new distribution obeys in a certain limit KNO scaling and that related n- 
particle correlation functions depend on two-particle correlation functions only. 
The explicit formula and its derivation are also given. Subsequent systematic 
analysis of final charged MD's in all available experiments in the accelerator 
region confirmed that NB MD is a very good candidate to describe MD in full 
phase spacal The main results of this accurate analysis have been presented 
at the Pavia Symposium in 1973. The analysis has been completed in 19711 
53 experiments were examined in total. The general trend of NB parameters 
is the same observed in cosmic ray physics, i.e., n and 1/k are growing as pi a b 
increases in all reactions involving protons on hadronic targets. In addition, 
at lower pi a b 1/k becomes negative after crossing the zero at approximately 30 
GeV/c. 

2.3 NB regularity in symmetric (pseudo) -rapidity intervals in pp energy do- 
main and then in all classes of reactions (Kyriat 1985, La Thuile 1989) 

NB regularity is rediscovered by UA5 collaboration and extended to symmetric 
pseudorapidity intervals in pp collisions at cm. energy ^/s — 200 GeV, 546 
GeV and 900 GeV. First results are presented by G. Ekspong at the Kyriat 
Symposium in 1985. After the impressive work by NA22, EMC, HRS and NA5 
collaborations - so well summarized by N. Schmitz in his talk at the La Thuile 
Symposium in 1989 - the new statement is: MD's in all classes of collisions 
in full phase space and in symmetric rapidity intervals are quite well fitted by 
NB MD. One problem shows up in pp at 900 GeV. NB does not fit data so 
well in the largest rapidity intervals. A shoulder is visible in the experimental 
distribution. 

2.4 Clan concept (Seewinkel 1987) 

Clan concept has been introduced in order to interpret the wide occurrence of 
approximate NB regularity in MP production. A clan is a group of particles of 
common ancestor (Sippe in German language). Each clan contains at least one 
particle (by definition). Clans are independently produced (by assumption). 
This new concept is discussed at the 1986 Seewinkel Symposium^ and exten- 



sively applied to experimental data in a series of papers by the same author!!. 
The distribution of particles in an average clan is logarithmic (or a weighted 
geometric) and its composition with independently emitted (i.e. Poissonian) 
clans leads to the observed NB MD. 



G NB (n,k;z) 



1 n(z-l) 



(3) 



G NB {n, k; z) = G Poisson [N; Gi og (n, k; z)] = exp [N (Gi og - 1)] (4) 

with N — fclog(l + n/k) and n c = n/N, the average number of clan and the 
average number of particles per clan respectively, i.e, Gjvs(n, fc; z) is a Com- 
pound Poisson Distribution (CPD) or discrete infinitely divisible distribution. 
At this level the introduction of the clan concept is essentially a change in the 
variables one suggests to observe in the multiplicity distribution, i.e. 

(:) * < 5 > 

Notice that 1/fc is interpreted here as a measure of aggregation of particles 
into clans 

i=m (6) 

i.e., it corresponds to the ratio of the probability to have two particles in 
the same clan, Pi (2), over the probability to have the two particles into two 
separate clans, ^2(2). 

Finally it should be recalled that in more general terms 1/fc is linked to the 
two particle rapidity correlations C 2 {y 1,2/2) and to the second order factorial 
cumulant, k 2 , by the equation 

-j; = K 2 = J C 2 (y 1 ,y 2 )dyidy2 (7) 

The deep connection among clan structure parameters and correlations is start- 
ing from these very simple remarks. 

The analysis of MD's in terms of the new parameters (clan structure anal- 
ysis) reveals new striking regularities: 

a. linear increase of ./V with the rapidity cut, y c , at fixed cm. energy; as 
y c approaches full phase space, N is bending towards a constant value; 

b. energy independence of N in a fixed rapidity interval in the region of 
linearity; 

c. iV is larger in e + e~ than in hh and in Ih reactions whereas n c is smaller 



2.5 NB regularity in Monte Carlo calculations (La Thuile Symposium 1989 
and then Perugia and Ringberg Workshops) 

In order to have informations on MD's at final parton level at different cm. 
energies in full phase space and in rapidity intervals , i.e., in a region where 
QCD predictions are poor, clan structure analysis has been applied to qq and 
gg systems as resulting in JETSET Monte Carlo calculations. Results of this 
search have been presented at the la Thuile Symposium in 1989 and in the 
same period at the Perugia and Ringberg Workshops. This study shows: 

a. final partons and charged hadron multiplicities in full phase space and 
in symmetric rapidity intervals are well fitted by NB MD's. 

b. Hadronic and partonic levels are not independent but linked by the 
following equations 



hadron 



parton 



= pn 

hadron parton 



(8) 



with p a constant <~ 2Qq/1 GeV, Qq is the parton virtuality cut-off, or, in 
terms of n-particle (parton) inclusive multiplicity distributions 



Qn 



(Vl, ■ ■■ ,Vn) = P n Qr 



hadron 



(yu- 



- Vn 



(9) 



parton 



which can be considered a generalization of preconfinement and defines the 
Generalized Local Parton Hadron Duality (GLHPD) hadronization prescrip- 
tion. 

c. The general trend of the average number of clans is at both levels 
consistent with what was previously discovered by clan structure analysis at 
hadron level. In addition the regularity is more impressive at parton level than 
at hadron level, i.e., where Altarelli-Parisi or Konishi-Ukawa-Veneziano QCD 
evolution equations predictions are not affected by the Monte Carlo hadroniza- 
tion prescription. 



2. 6 The shoulder structure puzzle and NB regularity in jets of fixed topology 
(La Thuile 1989, Momond 1991, Wuhan 1991, Vietri 1994, Stara Lesna 
1995) 

As already mentioned in pp reactions UA5 collaboration noticed that MD of the 
full sample of events at cm. energy y/s = 900 GeV in largest pseudorapidity 
intervals violates NB regularity and shows a characteristic shoulder structure. 
The effect was reported by Fuglesang at the La Thuile Symposium. Later the 
same effect was also seen at LEP at y/s — 91 GeV. 



The conclusion of those who were afraid of the regularity was that NB 
regularity was dead. 

At the same time more conservative people asked to look better at the 
data: shoulder effect might very well be due to the overlap of jets of different 
topologies i.e. the regularity should be tested at a more elementary level. 
The results of this new analysis in e + e~ at LEP energy are presented at the 
XXVI Rencontres de Moriond in 1991 by V. Uvarov. Although the analysis 
is dependent on the jet selection algorithm (JADE algorithm was used) three 
important facts occur: 

a. the shoulder in e + e~ annihilation is the result of the superposition of jets 
of different topologies; 

b. MD's in jets of different topologies are nicely fitted by NB MD's, i.e., 
the regularity which is violated in the full sample of events is approximately 
restored at a more elementary level of investigation; 

c. LEP data less approximate on sample of two jet events are fully consistent 
with HRS Collaboration data at lower cm. energy. 

In particular clan structure analysis reveals that standard trends for N 
and h c are less approximate on the sample of two jet events than in the full 
sample. In addition by applying GLHPD backwards one sees new impressive 
regularities at parton level: the average number of clan at parton level is 1 
in Ay =1, 2 in Ay=2, . . ., 5 in Ay=5. This fact is verified both at LEP 
and HRS cm. energies. The average number of clans is a linear function of 
rapidity variable. The average number of partons (particles) per clan follows 
the standard trend. It should be stressed that the same mechanism, i.e., the 
superposition of jets of different topologies (hard gluon radiation) explains in 
part also H q oscillations. 

2.7 New results (Faro 1996) 

Two new results have been presented at this Symposium, a. H q oscillations 
and NB regularity in quark jets of fixed flavour. In e + e~ annihilation at the Z a 
peak either the charged particle MD, P n , or H q oscillations are well reproduced 
in full phase space by the weighted superposition of two NB MD's associated 
to bb and light flavoured events, the weight being given by the fraction of bi 
events. The corresponding NB parameters are characterized by the fact thafl 

flight flavoured events 7^ ^bb events 5 ^ — (1^) 

Clan structure analysis suggests that the aggregation of particles into clan is 
the same in bb and light flavoured events! 



b. clan structure analysis at HI experiment at HERA. HI Collaboration at 
HERA presented data on MD's in different pseudorapidity windows at various 
energy intervals and then fitted these data by NB MEH. Clan structure analysis 
performed on these data reveals that the average number of clans in a fixed 
rapidity interval is energy independent and that the average number of particle 
per clan follows the same general trend already observed in EMC data (the 
figure presented at this Symposium is not shown here due to the lack of space) . 

3 Theory 

3.1 Jets as QCD Markov branching processes ( Goa 1979, Bruges 1980, Wuhan 
1991, Vietri 1994) 

QCD jets haveJ3een described in LLA as Markov branching processes of NB 
type since 1979u. The approach is a purely exclusive one and consists in solv- 
ing explicitly KUV equations with a fixed cut off regularization prescription, 
once the intrinsic Markoffian structure of the process has been recognized. It 
should be recalled that KUV equations are deduced in an inclusive approach to 
the production_arocess and are the differential formulation of the integral AP 
QCD equations!!! Results of this search have been extensively discussed in Goa 
(1979), Bruges (1980), Wuhan (1991) and Vietri (1994) Symposia. The main 
aspect of this study is that parton showers behave here as self-similar fractals; 
they develop according to a principle of disorder (fractal from Latin frangere, 
to break) followed by a principle of order and controlled by self-similarity, i.e., 
according to an "Italian firework" of jets of jets... In other words in a quark par- 
ton generated shower the ancestor quark "breaks" into bremsstrahlung gluon 
jets (the clans) which then decay into final parton via gluon self interaction, 
and partons belonging to each clan are distributed according to a weighted 
geometric distribution. What is interesting to point out is that clan structure 
analysis decouples the g — > g + g vertex from the q — > q + g vertex and that 
1/k is here just their ratio. 

3.2 Clan structure analysis in statistical physics (Aspen 1993, Vietri 1994) 
This study concerns: 

a. clans and compound Poisson distribution^; b. clans and hierarchical cor- 
relation functions^. It has been done in collaboration with Sergio Lupia and 
Roberto Ugoccioni and reported in Aspen 1993 and Vietri 1994 Symposia. 

a. The description of multiparticle production in a domain of rapidity Ay 
as a two step process can be formulated in very general terms in the frame- 
work of infinitely divisible distributions (IDD). Any discrete IDD is indeed a 



Compound Poisson Distribution (CPD). In the generating functions formalism 
this important property leads in a natural way to the generalized clan concept 
(g-clan) as can be seen by inspection of the following general definition of the 
generating function of any CPD: 

G C pd(z; Ay) = exp [N g (Ay)[G' (z; Ay) - 1]] (11) 

where N g is the average number of g-clans and G'(z;Ay) is the generating 
function of the MD inside an average g-clan. Since by definition each g-clan 
contains at least one particle the condition G' (z; Ay)\ z=0 = has to be im- 
posed. Accordingly clan concept is more general than NB multiplicity distribu- 
tion, whose wide occurrence in final particle multiplicity distributions justified 
its first introduction. In addition a simple relation exists for any CPD among 
the probability Po(Ay) to detect no particle in a given rapidity interval Ay 
and the average number of generalized clan, N g , i.e. 

Gcpd^ Ay) = e -^( A «) = P (Ay) (12) 

b. The main result can be summarized by the following theorem: factorial 
cumulant structure is hierarchical , i.e., n-order factorial cumulants, K n (Ay), 
are controlled by second order factorial cumulants, K2^Ay), only iff the function 
V(Ay) = — log Po ( Ay) /fi( Ay) scales with energy and Ay as a function of 
n(Ay)K 2 (Ay). Notice that the just mentioned scaling function is the inverse 
of the average number of particles per generalized clan, (fi c ) g . 

In conclusion clan structure parameters can be extended to the class of dis- 
crete IDD 's (they are the most economic way to describe a two step production 
process in which independently generated extended objects - clans, clusters, 
strings, ... - decay into final particles or partons). In addition these parame- 
ters are useful and intriguing in view of their deep physical meaning (they con- 
trol "voids" distribution in rapidity and being K n (Ay) — J A dy% . . . J Ay dy n 
c n (yi, ■ ■ ■ , Un) the hierarchical content of n-particle (-parton) correlation func- 
tions c n (yi, . . . , y n ) in the production process). Therefore the scheme described 
in (fsl) should be completed as follows 



\n c (Ay)J \V(Ay)J 



(13) 



In case of the NB MD the new set of parameters is linked to the previous one 
by the equations 



N(Ay) = -logPo(Ay) = - log 



k(Ay) 



n(Ay) + k(Ay) 



k(Ay) 



(14) 



The second equation sJapws that NB regularity predicts that n-particle corre- 
lations are hierarchical^. 

3.3 Clan structure parameters in theoretical models (Vietri 1994, Stara Lesna 
1995) 

The energy dependence of clan structure parameters can be calculated analyti- 
cally either in f.p.s. orJji rapidity intervals in the generalized simplified parton 
shower model (GSRS. 13. The model is a generalization of the work initiated 
with Leon Van Hovdi-l It is a two parameter model and is based on essentials 
of QCD (gluon self interaction) in a convenient kinematical framework: energy- 
momentum conservation laws are violated locally but not globally. Clans are 
considered here as independent intermediate fluctuating gluon sources. Ob- 
served general trends_pf clan structure parameter in energy and rapidity are 
correctly predictedElEj. These results suggest that the ancestors of the clans, 
i.e., the independently emitted bremsstrahlung gluons, should be indeed iden- 
tified at parton level with the just mentioned fluctuating gluon sources. 

3-4 Clans and the hadronization problem (Stara Lesna 1995, Faro 1996) 

At this point we are ready to ask the following question: are clans observable 
physical ohiects? A possible hint comes from the thermodynamical hadroniza- 
tion modebj. The model has been successfully applied to hadron production in 
e + e~ annihilation. The hadronic spectrum is quite well reproduced in terms of 
three parameters (the temperature, the volume and strangeness chemical sup- 
pression). The production process is described as a two-step process in which 
primary hadrons emitted by the thermal source decay into final particles. It 
is interesting to point out that final charged tracks MD's turn out to be of NB 
type and the average number of clans calculated from fitted NB coincides with 
the average number of primary hadrons predicted by the thermodynamical 
modeHl. This result suggests that at hadron level in e + e~ annihilation clans 
should be identified with primary hadrons. 

4 Conclusions 

The main results of experimental and theoretical work done along the years 
on approximate NB regularity and related clan structure, and reported in 
Multiparticle Dynamics Symposia, have been summarized. The hope is to 



have given a perspective for future work and a contribution to the answer of 
the following questions. 

Is NB regularity still an enigma - as Leon Van Hove pointed out in Shan- 
dong in 1987 - or a useful Arianna's thread - as I prefer now - in order to 
drive us in the labyrinth of multiparticle dynamics? 

Is the fact that in e + e~ annihilation NB MD can be associated to jets of 
different flavour accidental, or does it reveals the intimate structure of multi- 
particle dynamics in the building block of the reaction? If it is so, what about 
other reactions? 

Is clan concept a purely statistical one (like cluster expansion in statistical 
mechanics) or are clans observable physical objects? 

Are clans at parton level bremsstrahlung gluon jets? Can one consider 
the ancestors of each clan (the bremsstrahlung gluons) independently emitted 
fluctuating gluon sources? 

Are really clans at hadron level in e + e~ annihilation primary hadrons? If 
it is so, what about other reactions? 



Acknowledgements 

I thank Jorge Dias de Deus for the nice atmosphere created at the Conference. 
I also wish to express my gratitude to Sergio Lupia and Roberto Ugoccioni for 
a critical reading of the manuscript. 



References 



1. P.K. Mc Keown and A.W. Wolfendale, Proc. Phys. Soc. 89 (1966) 553. 

2. A. Giovannini, Nuovo Cim. 10A (1972) 713, ibid. 15A (1973) 543. 

3. A. Giovannini, Lettere al Nuovo Cimento 6 (1973) 514, A. Giovannini 
and M. Garetto ibid. 7 (1973) 35, A. Giovannini et al, Nuovo Cim. 24A 
(1974) 421. 

4. A. Giovannini et al., Nuovo Cimento 38A (1977) 38. 

5. L. Van Hove and A. Giovannini, in Proceedings of XVII International 
Symposium on Multiparticle Dynamics, Seewinkel, eds. M. Markitan et 
al., World Scientific, Singapore 1987, p. 561 

6. L. Van Hove and A. Giovannini, Acta Physica Polonica B19 (1988) 495, 
917, 931. 

7. A. Giovannini, S. Lupia an d R. Ugoccioni, preprint DFTT 50/96, MPI- 
PhT/96-80, LU TP 96-23, |hep-ph/9609306l (August 1996), to appear in 
Phys. Lett. B 

8. HI Coll. (S. Aid et al.), preprint DESY-96-160, subm. to Z. Phys. C. 



9. A. Giovannini, Nucl.Phys. B161 (1979) 429. 

10. K. Konishi, A. Ukawa, G. Vcncziano, Nucl. Phys. B157 (1979) 45. 

11. S. Lupia, A. Giovannini and R. Ugoccioni, Z. Phys. C 66 (1995) 195. 

12. S. Lupia, A. Giovannini and R. Ugoccioni, Z. Phys. C 59 (1993) 427. 

13. R. Ugoccioni, A. Giovannini and S. Lupia, Z. Phys. C64 (1994) 453. 

14. R. Ugoccioni and A. Giovannini, Z. Phys. C53 (1992) 239. 

15. A. Giovannini, S. Lupia, R. Ugoccioni, Z. Phys. C70 (1996) 291. 

16. F. Bccattini, Z. Phys. C69 (1996) 485. 

17. F. Becattini, A. Giovannini and S. Lupia, Z. Phys. C72 (1996) 491. 



